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Abstract 



For d > 1 and < (3 < a < 2, consider a family of pseudo differential operators {A Q + 
a /3 A /3/2. a G Qn R d tnat evo i ves continuously from A Q / 2 to A Q / 2 + A' 3 / 2 , ft gives arise 

C . to a family of Levy processes {X a , a <E [0, 1]} on R d , where each X a is the sum of independent 

a symmetric a-stable process and a symmetric /3-stable process with weight a. For any C • 
open set Del 11 , we establish explicit sharp two-sided estimates (uniform in a € [0, 1]) for the 
Qh transition density function of the subprocess X aD of X a killed upon leaving the open set D. 

The infinitesimal generator of X a - D is the non-local operator A Q + a^A^^ 2 with zero exterior 
condition on Z? c . As consequences of these sharp heat kernel estimates, we obtain uniform sharp 
Green function estimates for X a D and uniform boundary Harnack principle for X a in D with 
explicit decay rate. 
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1 Introduction 
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It is well-known that, for a second order elliptic differential operator C on M. d satisfying some 
natural conditions, there is a diffusion process X on M. d with C as its infinitesimal generator. The 
fundamental solution p(t,x,y) of dtu = Cu (also called the heat kernel of C) is the transition 
density function of X. Thus obtaining sharp two-sided estimates for p(t, x, y) is a fundamental 
problem in both analysis and probability theory. Such relationship is also true for a large class 
of Markov processes with discontinuous sample paths, which constitute an important family of 
stochastic processes in probability theory. They have been widely used in various applications. 
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One of the most important and most widely used family of Markov processes is the family 
of (rotationally) symmetric a-stable processes on M. d , < a < 2. A symmetric a-stable process 
X = {X t ,t>0,¥ x ,x £R d } on R d is a Levy process such that 



e i€-(Xt-x ) 



for every 16K and £ G 



When a = 2, X is a Brownian motion on M. d whose infinitesimal generator is the Laplacian A. 
When < a < 2, the infinitesimal generator of a symmetric a-stable process X in R rf is the 
fractional Laplacian A a / 2 , which is a prototype of nonlocal operators. The fractional Laplacian 
can be written in the form 

A"/y*) = A(d, -a) lim f (u(y) - u[x)) dy (1.1) 

£io 7{yead : |j/-x|>e} Ix-yl^ 

for some constant ^4(d, — a) := a2 Q ~ 1 7r~ d / 2 r(^±^)r(l — Here and in the sequel, we use := 

as a way of definition. Here T is the Gamma function defined by T(A) := f^°t x ~ 1 e~ t dt for every 
A > 0. 

Two-sided heat kernel estimates for diffusions in M. d have a long history and many beautiful 
results have been established. See (HUH] and the references therein. But, due to the complication 
near the boundary, two-sided estimates for the transition density functions of killed diffusions in 
a domain D (equivalently, the Dirichlet heat kernels) have been established only recently. See 
[12\ [T3l [T3] for upper bound estimates and [29] for the lower bound estimates of the Dirichlet heat 
kernels in bounded C > domains. In a recent paper [3], we succeeded in establishing sharp two- 
sided estimates for the heat kernel of the fractional Laplacian A a / 2 with zero exterior condition on 
D c (or equivalently, the transition density function of the killed a-stable process) in any C > open 
set. 

The approach developed in [3] can be adapted to establish heat kernel estimates of other jump 
processes in open subsets of M. d . In [4], the ideas of [3] were adapted to establish two-sided heat 
kernel estimates of censored stable processes in C 1 ' 1 open subsets of M, d . One of the main tools 
used in [3j is the boundary Harnack principle established in [2] and |17j . 

In [5] the ideas of [3] were adapted to establish two-sided heat kernel estimates of relativistic 
stable processes in C 1 ' 1 open subsets of R rf . One of main facts we used in [5] is that relativistic 
stable processes can be regarded as perturbations of symmetric stable processes in bounded open 
sets and therefore the Green functions of killed relativistic stable processes in bounded open sets 
are comparable to the Green functions of killed stable processes in the same open sets. 

The goal of this paper is to prove sharp two-sided estimates for the independent sum of an 
a-stable process and a /3-stable process, < (3 < a < 2, in C 1 ' 1 open subsets of M. d . Note that 
these processes can not be obtained from symmetric stable processes through a combination of 
Girsanov transform and Feynman-Kac transform. So the method of [5] can not be used to establish 
the comparability of the Green functions of these processes and the Green functions of symmetric 
stable processes in bounded open sets. Since the differences of the Levy measures of these processes 
and those of symmetric stable processes have infinite total mass, the method of [22J and [16J also 
could not be used to establish the comparability of the Green functions of these processes and the 
Green functions of symmetric stable processes in bounded open sets. The approach of this paper 
will be described in the second paragraph below after the statement of Corollary 11.21 
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Let us first recall some basic facts about the independent sum of stable processes and state our 
main result. 

Throughout the remainder of this paper, we assume that d > 1 and < (5 < a < 2. The 
Euclidean distance between x and y will be denoted as \x — y|. We will use B(x,r) to denote the 
open ball centered at x £ M. d with radius r > 

Suppose X is a symmetric a-stable process and Y is a symmetric /3-stable process on M. d and 
that X and K are independent. For any a > 0, we define X a by X t " := Xt + a^t- We will call the 
process X a the independent sum of the symmetric a-stable process X and the symmetric /3-stable 
process Y with weight a. The infinitesimal generator of X a is A Q//2 + a^A^/ 2 . Let p a (i, x, y) denote 
the transition density of X a (or equivalently the heat kernel of A a / 2 + a^A^/ 2 ) with respect to 
the Lebesgue measure on R d . We will use p(t, x, y) = p°(t, x, y) to denote the transition density of 
X = X° . Recently it is proven in |8j that 

At,x,y) * (r^Af^) A (^3^ + on (0, do) X M rf x R d . (1.2) 

Here and in the sequel, for a, b S R, a Ab := min{a, 6} and a V 6 := max{a, &}; for any two positive 
functions / and g, f X 5 means that there is a positive constant c > 1 so that c~ x g < f < eg on 
their common domain of definition. 

For every open subset D cR d , we denote by X a ' D the subprocess of X a killed upon leaving D. 
The infinitesimal generator of X a,D is (A Q / 2 + a@ A /3 ^ 2 )\e ) , the sum of two fractional Laplacians in 
D with zero exterior condition. It is known (see [8]) that X a,D has a Holder continuous transition 
density pp(t,x,y) with respect to the Lebesgue measure. 

Unlike the case of the symmetric a-stable process X := X°, X a does not have the stable scaling 
for a > 0. Instead, the following approximate scaling property is true and will be used several 
times in the rest of this paper: If {X^,t > 0} is the subprocess of X a killed upon leaving D, then 
{\- l X a x £,t > 0} is the subprocess of {xf (a ~ r))/l3 , t > 0} killed upon leaving A l D. So for any 
A > 0, we have 

pf}" D P)/ \t,x,y) = X d p a D (X a t,Xx,Xy) for t > and x, y € X' l D. (1.3) 
In particular, letting a = 1, A = a^/^ a ~^ and D = R d , we get 

Pd a/3 g /3 , 

p a (t, x, y) = at-Pp (a a -0t, a^i, a a -Py) for t > and x, y £ M • 

So we deduce from (|1.2|) that for any M > there exists a constants C > 1 depending only on 
d, a, (3 and M such that for any a £ (0, M] and (t, x, y) G (0, 00) x M d x M d 

C-V a (t, x, y) < p a (t, x, y) < C/ a (t, x, y), (1.4) 

where 

H«,.,») - ((^ A f -/«) A + ^^) . 

The purpose of this paper is to establish the following two-sided sharp estimates on p%(t, x, y) 
in Theorem 11.11 for every t > 0. To state this theorem, we first recall that an open set D in 
(when d > 2) is said to be a (uniform) C 1 ' open set if there exist a localization radius Rq > 
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and a constant Ao > such that for every z £ dD, there exist a C 1,1 -function 4> = cf> z : M.^ 1 — > R 
satisfying 0(0) = V0(O) = 0, ||V0||oo < Ao, |V0(x) — V(f)(z)\ < Aq\x — z\, and an orthonormal 
coordinate system CS Z with its origin at z such that 

B(z, R )nD = {y = (y, y d ) in CS Z : \y\ < R , y d > (j)(y)}. 

The pair (Ro,Aq) is called the characteristics of the C 1 ' 1 open set D. Note that a C ' open set 
D with characteristics (i?o 5 ^o) can be unbounded and disconnected; the distance between two 
distinct components of D is at least i?o- It is well known that any C 1 ' open set D satisfies both 
the uniform interior ball condition and the uniform exterior ball condition: there exists ro < Ro 
such that for every x € D with 5qd{x) < ro and y £M. d \D with Sgoiy) < ro, there are z^, z y E <9.D 
so that |x — z x \ = 5g£>(x), \y — z y \ = Sgrjiy) and that B(xo,ro) C D and B(yo,ro) C M. d \ D for 
= z x + ro(a; — — z x \ and yo = % + r o(y — — %|- By a C open set in R we mean 

an open set which can be written as the union of disjoint intervals so that the minimum of the 
lengths of all these intervals is positive and the minimum of the distances between these intervals 
is a positive constant Rq. 

Theorem 1.1 Suppose M > 0. Let D be a C 1,1 open subset ofW d and Sjj(x) the Euclidean distance 
between x and D c . 

(i) For every T > 0, there is a positive constant C\ = C±(D, M,a, (3,T) > 1 such that for every 
a e (0,M], 

Ci" 1 mt,x,y) <p a D {t,x,y) <dmt,x,y), 

where 




Vi J \ Vt 

(ii) Suppose in addition that D is bounded. For every T > 0, there is a constant Ci = C2(D, M, a, (5, T) > 
1 so that for every a S (0, M] and (t, x, y) E [T, oo) x D x D, 

C^ 1 e-^ 1 5 D {xY' 2 5 D (y) a l 2 < P a D (t,x,y) < C 2 e" Alt 5 D {x) a l 2 5 D (y) a / 2 , 

where X\ > is the smallest eigenvalue of — (A Q//2 + a /3 A@/ 2 )\rj- 

Letting a — > 0, Theorem 11.11 recovers the heat kernel estimates for symmetric a-stable processes 
obtained in [3]. By integrating the two-sided heat kernel estimates in Theorem 11.11 with respect 
to t, we obtain the following estimates on the Green function G a D {x,y) := J °° pp(t,x,y)dt, which 
mean that, for bounded C ' domains D, G a D and G® D are comparable, see [Sj and [20J. To the best 
of our knowledge, the Green function estimates in the corollary below are new. 

Corollary 1.2 Suppose M > 0. For any bounded C l,1 -open set D in M. d , there is a constant 
C 3 = C 3 (D, M, a,P)>l so that for every a £ (0, M], 

CjT 1 9D(v,y) < G a D (x,y) < C 3 g D {x,y) for x,y £ D, 
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where 

' (i A SD{xr ^ yr/2 ) \x - y\ a - d when d > a, 

9D{x,y) := I log (l + SDix) "Z 5 y^ y)a/2 ) whend = l = a, (1.5) 

(5 D {x)5 D {y)) K A " v ; — whend = l<a. 

Theorem ll.l( i) will be established through Theorems 12.81 and 13.51 which give the upper bound 
and lower bound estimates, respectively. Theorem 11,1 H i) is a consequence of the intrinsic ultracon- 
tractivity of X a in a bounded open set and the continuity of eigenvalues proved in [10J. In fact, 
the upper bound estimates in both Theorem 1 1 . 1 1 and Corollary 11.21 hold for any open set D with a 
weak version of the uniform exterior ball condition in place of the C ' condition, while the lower 
bound estimates in both Theorem 11.11 and Corollary 11.21 hold for any open set D with the uniform 
interior ball condition in place of the C ' condition (see Theorems 12.81 and 13.51 and the proofs for 
Theorem II .lf fi) and Corollary 1 1 . 2[) . 

Although we follow the main ideas we developed in [3], there are several new difficulties in 
obtaining two-sided Dirichlet heat kernel estimates for X a : Even though the boundary Harnack 
principle has been extended in [19] to a large class of pure jump Levy processes including X a , the 
explicit decay rate of harmonic functions of X a near the boundary of D was unknown. Instead, 
following the approach in [6], we establish necessary estimates using suitably chosen subharmonic 
and superharmonic functions of the process X a . As in [6], we need to use finite range (or truncated) 
symmetric /3-stable process Y x obtained from Y by suppressing all its jumps of size larger than A. 
The infinitesimal generator of Y x is 

m\{x) := A(d, -a) lim / (u(y) - u(x)) dy (1.6) 

J { y m d :e<\y-x\<\} \ x ~ V\ H 

When A = 1, we will simply denote A^ 2 by A@/ 2 . We first establish the desired estimates for the 
Levy process X a := X + aY x l a . The infinitesimal generator of X a is A°/ 2 + a^A^/ 2 . The desired 
estimates for X a = X + aY can then be obtained by adding back those jumps of Y of size larger 
than 1/a. To obtain the lower bound of p a (t,x,y), we use the Dirichlet heat kernel estimate for 
the fractional Laplacian in [3] and a comparison of the killed subordinate stable process with the 
subordinate killed stable process where we will use some of the results obtained in [26] . 

We like to point out that unlike [3] the boundary Harnack principle for X a is not used in this 
paper, which indicates that it might be possible to obtain sharp heat kernel estimate for processes 
for which the boundary Harnack principle fails. 

As a consequence of Corollary 11.21 we have the following uniform boundary Harnack principle 
with explicit decay rate. 

Theorem 1.3 Suppose that M > 0. For any C 1 ' 1 open set D in M. d with the characteristics (Rq, A), 
there exists a positive constant C4 = C^a, /3, d, A, R ,M) > 1 such that for a E [0, M], r E (0, Rq], 
Q E dD and any nonnegative function u in W d that is harmonic in D n B(Q,r) with respect to X a 
and vanishes continuously on D c n B(Q,r), we have 

44<C4^^ for every x,y e D (1 B(Q,r/2). (1.7) 
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Throughout this paper, we will use capital letters C\, C2, . . . to denote constants in the state- 
ments of results, and their labeling will be fixed. The lower case constants ci,P2,... will denote 
generic constants used in proofs, whose exact values are not important and can change from one 
appearance to another. The labeling of the lower case constants starts anew in every proof. The 
dependence of the constants on dimension d may not be mentioned explicitly. For every function 
/, let / + := / V 0. We will use d to denote a cemetery point and for every function /, we extend 
its definition to d by setting f(d) = 0. We will use dx to denote the Lebesgue measure in R rf . For 
a Borel set A C M. d , we also use |A| to denote its Lebesgue measure. 

2 Upper bound estimate 

Throughout this section we assume that D is an open set satisfying the uniform exterior ball 
condition with radius ro > in the following sense: for every z £ dD and r £ (0, ro), there is a ball 
B z of radius r such that B z cM d \D and dB z n dD = {z}. The goal of this section is to establish 
the upper bound for the transition density (heat kernel) Pjj(t,x,y). One of the main difficulties of 
getting the upper bound for p a D {t,x,y) is to obtaining the correct boundary decay rate. 

Recall that A a / 2 and A^ 2 are defined by (jl.ip and (|1.6|) . The next two lemmas can be proved 
by direct computation, whose proofs can be found in [17] and [6], respectively. 

For p > 0, let w p (x) := (x~l) p . 

Lemma 2.1 For any x 6 (0, 00) x R* -1 , we have 

A a/2 w a/2 (x) =0. (2.1) 

Moreover, for every p £ (a/2, a), there is a positive constant C5 = C$(d, a,p) such that for every 
x e (0,oo) x R d ~ l 

A a / 2 w p (x) =C 5 x?" a . (2.2) 

Lemma 2.2 There are constants i?* £ (0, 1), Cq > Cj > depending onp, d and a only such that 
for every x £ (0, i?*] x IR^ -1 

C 7 x\' a < A a / 2 w p (x) < CQx\~ a for a/2<p< a, (2.3) 

\A a / 2 w p (x)\ < C 6 I log xi I forp = a (2.4) 

and 

\A a / 2 w p (x)\ < C 6 forp>a. (2.5) 

In the remainder of this paper, R* will always stand for the constant in Lemma 12.21 The 
following result and its proof are similar to Lemma 3.2 of [6] and the proof there. For reader's 
convenience, we spell out the details of the proof here. 

Lemma 2.3 Assumethatr 1 G (0,1/2] andp > f . Let Si := i?*A(n/4), U := {z £ R d : n < \z\ < 3n 
and 

h P (y) ■■= (yd - \jr\ - \y\ 2 ^ lun{z d >o,\z\< ri /2}{y)- 
Then there exist Ci = Ci(a,p, r\) > 0, i = 8, • • ■ , 12, such that 
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(i) when p G (a/2, a), we have for all y G {z G M d : z d > 0, ri < |z| < ri + Si, \z\ < ri/4}, 

c 8 L - V^-lyl 2 ) P ° < A a/2 My) < c 9 (yd - V r i-lyl 2 ) P ° ( 2 - 6 ) 

and 

C 8 - ^ - |y| 2 ) P < A^ 2 h p (y) < C 9 (y d - yj r\ - |y| 2 ) " ° ; (2.7) 

(ii) w/ien p > a, we /iai>e 

|A Q/2 /i p (y)| < Cio /or a// y £ {z € R d : z d > 0, n < \z\ <r x + 6 1: \z\ < n/ij ; (2.8) 

(iii) when p = a/2, we have 

\A a / 2 h a/2 (y)\ < Cii /or all y e {z e R d : z d > 0, n < \z\ <n + S, \z\ < n/4} ; (2.9) 

(iv) when p = a, we have for every y G jz G R d : z^ > 0, r\ < \z\ < n + 5, |2j < ri/4}, 

|A Q / 2 v 2 (y)l < Ci2 

Proof. Let 

r(y) := -\/r 2 - |y|2 and h(y) := y d - T(y), y G U. 

Fix x G {z G M d : z d > 0, ri < |z| < ri + i?* A (ri/8), |i] < ri/4} and choose a point xo G 
9i?+(0, n) := {zrf > 0, \z\ = 7*1} satisfying x = 2o. Denote by n (xq) the inward unit normal vector 
at xq for the exterior ball B(0,r\) c and set <J?(y) = (y — XQ, ~n(xo)) for y G M d . IT = {y : <&(y) = 0} is 
the plane tangent to dB + (0, r\) at the point x$. Let r* : x G M^ -1 — > R be the function describing 
the hyperplane IT, that is, ((x,T*(x)) — xq, Ti(xq)) = 0. We also let 

E := {y = (y, y d ) ■ y G U, \y - x\ < ri/4} , 
A := {y : r*(y) > y d > T(y), |y - x\ < n/4} 

and h(y) := (y d - T*(y))l^ yd>v *^(y) for y G M d . Since VT(x) - Vr*(x) = 0, by the mean value 
theorem 

1%) - Hy)\ < |r(50 - r*(y)| < A|y - x\ 2 , y eE. (2.11) 

Let £ n (y) = dist(y,II) for y G R d and £7 r * = {y G R d : y d > T*(y)}. Let 6;, := y 7 ! + \VT(x)\ 2 and 

h x , p (y) := (h(y)T. 

Note that h x , P (x) = h p {x) and B(x,r\/4) C\TJ <Z E. Since /i(y) = b x 5 n (y) on Z?r*, by Lemma [23] 

A a/2 ^, a/2 (x) = (2.12) 
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and, if a/2 < p < a, 



A a / 2 h XiP (x) = Cl bPJP- a (x) = Cl b a x {h{x)y- a 



(2.13) 



for some c\ > 0. By Lemma \2~2[ there are constants q > 0, z = 2 ... 6, such that for y € Dp* and 
5 n (y) < i?*, when a/2 < p < a, 

c 2 (h(x)f- a < c 3 FJ p n - a (x) < A a / 2 h x>p (x) = bP x A a / 2 (5 n (x)r < c 4 ^- a (x) < c 5 (h(x)f- a , 

(2.14) 



when p > a, 



and when p = a, 



Note that 

\A a / 2 (h p -h x , p )(x)\ = A{d,-a) 

< A(d, -a) 



\A a / 2 h x , p (x)\ = V x \A a/2 {5 Il {x)Y\ < c 6 
\A a ' 2 h x ^x)\ <c 6 \log(h(x))\. 

(h p {y) - h P)X {y)) 



(2.15) 
(2.16) 



lim 



dy 



+ A(d, —a) lim 



40 J{i>\ y -x\>s} \ x - y\ d+a 

(h p (y) - h p , x {y)) 

i i d+a y 

{l>\y-x\>n/4} \ x ~ U\ 

\h P (y) - h P)X (y)\ 



(2.17) 



< cj + A(d, —a 
= : c 7 + h+ h 

and, similarly, 

\A a / 2 (hp-h x , p )(x)\ 
<A(d, -a 



10 J{ ri /4>\y-x\>e} \% ~ V\ d 

h P (y) + h P)X (y) 



dy 



\x - y\ 



d+a 



dy + A(d, —a) 



\h P {y) - h PtX (y)\ 



\ x - y\ 



d+a 



(2.18) 



+ A(d, —a 



{\y-x\>r±/4} \ x V\ 

h P (y) + h P;X (y) 



(My) - hAy)) d 
d+a y 



dy + A(d, —a) 



\h P (y) ~ h pAy)\ 

E \x-y\ d + a 



\x - y\ d+a 
Since for y 6 B(x, ri/4) c , 

\h x , P {y) - K,p( x )\ < c 8 \y - x\ p and \h p (y)\ < c 8 
and h p (y) = for \y\ > r\/2, for a/2 < p < a we get 



:h + h+h- (2.19) 



I z <A(d, -a 
+ A(d,-a) 

<C9 / 

JB(x,ri/A) c 



\h x , P (y) - h x , p (x)\ 



B(x,n/4)° \ x ~ y\ d+a 



dy + A(d, —a) 



\h P (y) - h p (x)\ 



S(x,ri/4) c n{|y|>ri/2} \ x ~ 2/1 



h p {x) 

—iTa- d y 



1 



\x - y\ 



d+a—p 



-dy + c 9 



B(x,r-,_/4) c 



\ x - y\ 



d+a 



B(x,n/4:) c r){\y\<r 1 /2} \ x ~ 2/1 



dy < cio < oo. 



d+a 



dy 



(2.20) 
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We claim that, if p > a/2, 

h + h < en < oo. (2.21) 

Note that for y £ A 

\kcM\ + \h p (y)\ < \y d - r*(y)\*> + \y d - T(y)f < 2\T(y) - T*(W 

< 2\T(y) - F(x) - VT(x) ■ (y-x)\ p < 2^ 12 \y-x\ 2p . (2.22) 

Furthermore, since |T(y) — T*(y)| < c\^\y — x\ 2 < C\2T 2 on \y — x\ = r, this together with (|2.22p 
yields that 

h < c u [ ri/ \ 2p - a - d [ l A (y) m d - x {dy)dr 

Jo J\y-x\=r 

' l/ \ 2p - a - d m d ^{{y. \y-x\=r, V* (y) > y d > T(y)})dr 



Cl4 





rn/4 



< c i5 / r 2p ~ a dr < oo. 
Jo 

Note that for y £ E 

\h P (y) - h XjP (y)\ <c 16 \(h(y)r - (h(y)Y\ < c 17 (h(y))^-\h(y)-h(y)\, (2.23) 
where (p — 1)_ := {p — 1) A 0. In the last inequality above, we have used the inequalities 

\V -aP\ < V } - l \b- a\ for a,b > 0, <p < 1 

and 

\V ~ aP\ <(p + l)\b- a\ for a, b G (0, 1), p > 1. 

For y = (y,y d ) € M. d , we use an affine coordinate system z = (z,z d ) to represent it so that 
z d = Ud~ r* (y ) and z is the coordinates in an orthogonal coordinate system centered at xq for the 
(d — l)-dimensional hyperplane II for the point (y, T*(y)). Denote such an affine transformation 
y i— ► z by z = ^(y). It is clear that there is a constant cis > 1 so that for every y £ M. d , 

C 18 \V -%\ < \z\ < ci 8 \y-x\, c^\y - x\ < - < c 18 \y - x\ 

and that 

^(E) C {z = (z, z d ) G R d : \z\ < cisn and < z d < Cisn}. 

Denote x d — T*(x) by w; that is, ^(x) = (0,w). Hence by (|2.1ip and (|2.23j) and applying the 
transform we have by using polar coordinates for z on the hyperplane II, 

h < ci9 / -^n r±- dy < c 19 / — 2__ — !_L — dz 



ci$ri / /-cisri r d-2 

(P-1)- 







< C20 / \ t 1 tttt: T<^ r 

/- d Wo (r + |z d -u;|) d + Q - 2 / 







Cl8ri (p-1) ( 1 1 



< C 2 2 / 773 dz d < C 23 < OO, 

'o z\ p> \z d -w\ a - 1 



9 



where all the constants depend on a, p and r\. The last inequality is due to the fact that since p > 0, 
< a < 2 and (1 — p) + + a — 1 = max{a — p, a — 1} < 1, by the dominated convergence theorem, 



w 



) == Jo 



Cl8^1 



1 



hence is bounded. 



dzd is a strictly positive continuous function in Xd G [0, cisri] and 



Thus we have proved the claim (|2.2ip . The desired estimates (|2.6[) - (|2.10j) now follow from 



(I2T21 - (|2T2T1) . 

It is well-known that X 1 has Levy intensity 



□ 



JH*,v) = jH\*-v\) = ^^+ A{d '-® 



A scaling argument yields that 



J a (x,y)=f(\x-y\) 



x _ y]d+a ^ x _ 



Aid, -a) , A(d, -P) 
~\~ ci 



\x - y\ 



\x-y 



\d+f3- 



Put 



(2.24) 



r(r):=l + aP A } d ; S- r a ^, r > 0. 
A(d, —a) 

Clearly for a G (0, M] and r > 0, 1 < ^ a (r) < 1 + cM^r a ~ p and 

J a (x, y) = f(\x -y\) = A[dl ~ d a } a V{\x - y\). 

\X — y\ a + a 

The Levy intensity gives rise to a Levy system for X a , which describes the jumps of the process 
X a : for any non- negative measurable function / on M + x M. d x M. d , x G M. d and stopping time T 
(with respect to the filtration of X a ), 



s<T 



E, 



T 



v JR d 



f(s,X^y)J a (X^y)dy)ds 



(2.25) 



(See, for example, Proof of Lemma 4.7] and [SJ Appendix A]). 

For any open set D C TSL , let = T a (D) := inf{i > : X" ^ D} denote the first exit time 
from D by X a . 

The next lemma follows immediately from a special case of |19t Proposition 2.10 and Lemma 
3.6] and the scaling property of X a . 

Lemma 2.4 For any 6, M G (0,oo), there exists C13 = C13 (M, b, a, (3) > such that for every 
x G M. d , a G [0, M] and r G (0, 6], 



B(x ,r) 



< C13 r a/2 (r - |s-so| 



\a/2 



/or x G B(xq, r). 



(2.26) 



For A > 0, y A = (l^\PxO is a Levy process in R such that 

Ex \ e W x -y x )] = e ~tm for every x G M d and £ G M d , 
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with 



i>(0=A(d,-(3) 



1 - cos(£ • y) 



\d+p 



dy. 



'{|J/|<A} 

In other words, Y x is a pure jump symmetric Levy process in M. d with a Levy density given by 
A(d, — (3)\x\~ d ~P 1{u|<a}- For a > 0, suppose Y l / a is independent of the symmetric a-stable process 
X on M. d . Define 



X?:=X t + aY t 



1/a 



t > 0. 



We will call the process X a the independent sum of the symmetric a-stable process X and the 
truncated symmetric /3-stable process Y x l a with weight a > 0. The infinitesimal generator of X a 
is A"/ 2 + a^l 2 . 

For any open set U C M d , let = inf{t > : Xf ^ U} be the first exit time from U by X a . 
The truncated process X a will be used in the proof of next lemma. 

Lemma 2.5 Assume r\ G (0, \] and M > 0. Let t7 := {z G M d : n < |z| < 3ri/2}. T/iere are 
constants C14 = Ci4(ri, a) > and C15 = Cis(ri, M, a, (3) > suc/i /or every a £ [0, M] 



E 



K] < Ci 4 P. > 3n/2) < Cis^^)"/ 2 , /orn < |z| < 5 n /4. 



(2.27) 



Proof. The first inequality in (|2.27p is easy. In fact, by the Levy system (|2.25j) with 

f(s,x,y) = l[/(x)l {5ri< |j / | <10ri} (y) 



and T = Ty, we have that for x G U 



X? a \ > 3n/2 > 10ri > \X$,\ > 5n 



> E, 



J{5ri<|y|<10ri} 



r(X^,y)dyds 



«4(cZ, —a) 
{5n<|»|<ion} -y| d+a 



> ciE x [t£], 



where ci = c(ri, a) > 0. 

It is enough to prove the second inequality in (|2.27p for n < \x\ < n + 5 for some small 5 > 0. 
Without loss of generality, we assume x = and x^ > 0. Let p > be such that p ^ (5 and 

a - (/3/2) < p < (a - (0/2) + (a - /3)/3) A a. 

Note that a/2 < p < 3a/2 - (5. Define 



2/d 



a/2 



\l r \- W) lun{z d >0,\z\<n/2}{y), 



9 P (y) ■= (yd - \Jr\ - \y\A lun{z d >o,\z\<n/2}(y), 



and let <fi be a smooth function on M. d with bounded first and second partial derivatives such that 
0( y ) = 2 4+ P\y\ 2 /rj for y G {z rf > 0,n < \y\ < 4n/5, |z| < n/4} and 2? < < 4? if |y| > n/2 
or |y| > 3ri/2. 
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Since ri < 1/4, it is easy to see that ||g p ||oo < 1- Now we define 

u(y) ■= Ky) + 4>(y) - 9 P (y)- 

By Taylor's expansion with the remainder of order 2, we get that for any a G (0, M] and y G M. d , 
|(A a/2 + a /3 A /3/2 )</>(y)| < || ^" /2 ^|| 00 + M p \\ A^^L < c 2 (a,/3,M) < oo. (2.28) 

Moreover, by (p76]) - (pT8|) . there exist c 3 = c 3 (a,/3) > and 5\ = 8\{a,0) G (0,ri/8) such that 
& a/2 g P (y) > C3S D (y) p - a for ye{zeR d :z d >0, ri < \z\ < r x + 5 U \z\ < n/4} 

and 

A P/2 9 P (y) > -c 3 5 D (y)( p -V A0 for j, G {z G K d : z d > 0, n < \z\ < n + <5 X , \z\ < n/4} . 

Thus there exist c 4 = 04(0, (3,M) > and 52 = $2 (a, /?, M) G (0, 61) such that for all a G (0, M] 
and y e {z eR d : z d > 0, ri < \z\ < T\ + <5 2 , |i] < n/4}, 

(A«/ 2 + a^l 2 )g p {y) > c 3 5 D ( y y- a - C3 M^ D (y) {p ^ )A0 > c 4 5 D ( y y- a . (2.29) 

Furthermore by §2J)§ and (f2jg jnf210|l . there exist c 5 = c 5 (a,/3,M) > and 5 3 = 5 3 (a,/3) G (0, Si) 
such that for all a G (0, M] and for every y G {2: G M rf : z d > 0, n < \z\ < T\ + #3, |J| < n/4} , 

(A Q / 2 + a A^ 2 )h(y) < \ A a ' 2 h{y)\ + M p \ A^ 2 h(y)\ 

< fcs + csfe^)^ 2 "^ if a/2, 
\c 5 + c 5 \logS D (y)\ if /3 = a/2. 

Since p - a < a/2 - /3, by ([2^ - (j2~29j) . there exists 64 = <5 4 (a, /?, M) G (0, <5 2 A <f 3 ) such that for 
all a G (0,M] and y e V := {z e R d : z d > 0, n < \z\ < n + &t, |S] < n/4} 

(A«/ 2 + aPW 2 )u(y) < c 2 + c 5 + c 5 («M2/) (a/2_/3)AO + | log Ml/) f) - c 4 My) P_a < 0. (2.31) 



Let t] be a non-negative smooth radial function with compact support in R rf such that rj(x) = 
for |x| > 1 and J" Kd rj(x)dx = 1. For fe > 1, define rjk(x) = 2 kd rj(2 k x). Set u^ k \z) := (% * u)(z). As 
(A Q / 2 + a^A 3 / 2 )^ =r] k * (A a / 2 + a /3 A /3 / 2 )n, we have by (I23TT) that 

on V k := {z e R d : z d > 0, n + 2~ fe < |z| < n + £4 — 2~ fc and |£| < n/4 - 2~ k }. Since u( fe ) is a 
bounded smooth function on M. d with bounded first and second partial derivatives, by Ito's formula 
and the Levy system f|2.25[) . 

M k := u( fc )(A7) - n( fc )(X a ) - J* (A a / 2 + (/A 3 / 2 ) u«(X s a )ds (2.32) 

is a martingale. Thus it follows from (|2.32|) that t \—* (X?" A ~ a ) is a bounded supermartingale. 
Since V k increases to V and u is bounded and continuous on V, we conclude that 



1 1— ► uiXl^a I is a bounded supermartingale. (2.33) 
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We observe that, since (j)(x) = 0, 

u(x) < 8u(x) a/2 . 

We also observe that, since <f> > 2g p outside of {z £ U : Zd > 0, \z\ < d/2} and 

a/2 



(2.34) 



u{y) > (yd - \f\ 



yd 



> c 6 



on {zd > 0, n + 04 < | | < 3ri/2, \z\ < d/2}, we have 



(2.35) 



u{y)>c 7 >0 for y £V c \B(0, ri ), 

where cj depends on o~4, a, (3 and r%. Therefore, by (12,33j) - (|2.35j) we get 

5u(x) a/2 > u(x) > E x . > C7 P X . (x$a € V c \ 5(0, n)) > c 7 F x (\X^\ > 3n/2) . (2.36) 

Note that there exist c$ = c$(a, d, d) > and eg = cg(/3, d, ri) > such that for z £ U, 

dy ^ _ f dy 



< c 8 



^2 



and 



{|y|>2n} \ Z ~ V\ d+a J{2r 1 <\y\<3r 1 } \ z ~ V\ d+C 

dy ^ f dy 



< Cg 



'{|»|>2n} \ z ~ y\ d+P J{2rx<\y\<3n} \ z ~ y\ d+f3 ' 
Thus by (|2.25|) . there exists a positive constant do = c\o(d, a, (3, M) such that for any a € (0, M], 



\X%\ >2n = 



J{\y\>2n} 



J a {X a s ,y)dyds 



10 J{2n<|j/|<3ri} 

ci P a! (3ri>|^|>2n). 



J a (X*,y)dyds 



(2.37) 



Since r% < 1/4 and the processes X and Y do not jump simultaneously, we have by (|2.36|) that 
there is a positive constant di = cn(d,a, (3, M,r±) such that for all a € (0, M], 

F x (\X?*\ > 3n/2) < (do + 1)1?. (3d > |X«a| > 3d/2) 
= (do + (3d > > 3d/2) 

< (c w + l)F x (\X$a\ > 3d/2) < c u Su(x) a/2 . 

□ 

Lemma 2.6 Assume M > and d G (0, |]. Ze£ £ = {iel d : |x| > d}. T/ien /or every T > 0, 
i/iere is a constant C±q = C%q (t\ , a, /3, T, M) > suc/i t/iai for every a E [0, M], 



p%(t, x, y) < Ci 6 ^(x) Q/2 J a (x, y) /or n < |s| < 5n/4, \y\ > 2 n and t < T. 
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Proof. Define U := {z 6 R d : n < \z\ < 3ri/2}. Since X a satisfies the hypothesis H in [28], by 
[281 Theorem 1], X"a dU with probability 1. For r x < \x\ < 5ri/4, \y\ > 2r\ and t G (0,T], it 
follows from the strong Markov property of X a and (|2.25p that 

p%(t,x,y) = E x [p a E (t-T%,X? tr ,y) : r% < t 
= I / pu(s,x,z) I / J a (z, w) p%(t — s, w, y)dw \ dzds 

Jo JU \J{w:\w\>3r 1 /2} J 

= I / pu(s,x,z)\ / J a (z,w)p%(t — s,w,y)dw \ dzds 

Jo Ju \J{w:(3r 1 /4)+(\y\/2)>\w\>3r 1 /2} ' J 

+ / pu(s,x,z)\ / J a {z, w)p a E {t — s, w, y)dw \ dzds 

Jo Ju \/{iu:|iu|>(3ri/4)+(| W |/2)} " J 

=:I + II. 

Note that for \w\ < (3ri/4) + (|y|/2), 

I»- 9 |>| S |-H>5(M-^)>|1>^. (2.38) 

Since p%(t — s,w,y) < p a (t — s,w,y), by (|1.4[) and (|2.38p . there exist constants c\ = ci(a, /3, M) > 
and c 2 = c 2 (a, /3, M) > such that for a G (0, M] 

I< f f p1f{s,x,z)l [ J a (z,w)c 1 TJ a (w,y)dw) dzds 

Jo Ju \^{m:(3n/4)+(|v|/2)>|iu|>3ri/2} / 

<c 2 TJ a (x,y) [ [ pfj(s,x,z)l I J a (z,w)dw) dzds 

Jo Ju W{w:3|x-2/|/4>|w|>3n/2} / 



= c 2 TJ a (x,y)¥ x (3n/2 < < 3\x-y\/4; r% <t] 

<c 2 TJ a (x,y)¥ x (\X%\ >3ri/2). 

By Lemma 12.51 . we have for |x| G (ri,5ri/4), 

> 3n/2) < c 3 5u(x) a / 2 = c 3 S E (x) a / 2 

for some positive constant C3 = cs(r±,a,/3,M). Thus 

/< c 4 (TVl)^(x) Q/2 J a (x,y) (2.39) 

for some positive constant C4 = c^(r\,a,, j3,M). On the other hand, for z G £/ and ty G M rf with 
M > (3n/4) + (|y|/2), 

I* - -I > M - N > - ( \y\ - ^ V M > tM. 

1 1 - 1 1 2 V 2 - 8 - 16 
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Thus by the symmetry of p%(t — s, w, y) in (w,y), we have 

II <c 5 J a {x,y) [ I pfj (s, x,z) ( / p%(t- s,y,w)dw ) dzds 

JO JU \ J (w:\w\X3n /4)+(\y\/2)\ J 



<c 5 J a {x,y) / / p u (s,x,z)dzds 
Jo JU 



= c 5 J a (x,y)E x [T u ] < c G 5 E (x) a / 2 J a (x,y) 

for some positive constants Ck = Ck(ri,a,(3,M), k = 5, 6. In the last inequality, we used Lemma 
12,51 to deduce that K x [tjj] < c5u(x) a / 2 = c<5^(x) a//2 for some positive constant c = c(r±, a, (3, M). 
This together with (12. 39ft proves the lemma. □ 



Theorem 2.7 Assume that M > and D is an open set that satisfies the uniform exterior ball 
condition with radius ro > 0. Then for every T > 0, there is a constant Cyj = Cn(ro/T, a, f3, M) > 
such that for all a G (0, M], A G (0,T] and x,y £ \~ l D, 



a/2 



p a x _ lD (l,x,y) < Ci7 (1 A J a (x, y)) 6 x -i D (x) 



Proof. Note that for every A G (0, T], A 1 D satisfies the uniform exterior ball condition with 
radius r /T. For x,y G A" 1 !), let z G 5(A _1 Z?) be that \x - z\ = S x -i D (x). Let B z C (A _1 L>) C be 
the ball with radius r x := 4" 1 A (r Q /T) so that &B 2 n 9(A _1 L>) = {z}. Since, by ([13]) 

p a x ^ lD {l,x,y) <p a {l,x,y) <c(lAJ a {x,y)), 

it suffices to prove the theorem for x G \~ X D with < ri/4. When (Ja- 1 !)^) < r i/4 and 

| a; — 2/| > 5ri, we have > 2ri and so by Lemma [231 there is a constant c\ > that depends 

only on (ro/T, d, a, /?, M) such that for t G (0, 1], 

p$-i D (t,x,y) < pffe M)c (t,x,y) < c 1 5 (Ez)c {x) a / 2 J a {x,y)=c l 5 x -, D {x) a l 2 J a {x,y). (2.40) 

So it remains to show that when o~\-\d{x) < r i/4 and |x — y| < 5ri, there exists a positive constant 
C2 = C2(ro/T, d, a, (3, M) such that 



(2.41) 



Let z x G <9(A 1 D) be such that \x — z x \ = o~\-id{ x ) an d z o £ K so that 

5(zo, n) C (A~ 1 L») C and dB(z , ri ) n ^(A" 1 ^) = {z^}. 
Define [/:={wel d : |w - z | G (n, 8ri)}. Note that 

x,y E U H \~ 1 D and <5;/(:e) = 5a-i_d(x)- 
By the strong Markov property and the symmetry of p x _ i D (l, x, y) in x and y, we have 



p a x - 1D (l,x,y) = p a UnX ^ D (l,x,y)+E y p a x -i D (l - r% nX - lD , X% 



< 1 
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By the semigroup property and JL 



J&nA-^CM.J/) = / Pun\-^D( l / 2 i x i z )Pun\-^D( l l 2 ' z iy) dz 
Jun\- 1 D 

< ||p ffl (l/2,-,.)||ooP4^nA-iD>l/2) 

< C 4 ^(x) a/2 = C 4 <5 A - 1D (x) a / 2 . 

In the last inequality, we used Lemma 12.51 

On the other hand, we have li G U c n A -1 !) on {i\ a , . lf] < 1}, and so 



|-X£> - X \ 

un\- 1 D 



> 7n, on {r^ nA _ 1D < 1} . 
Consequently, by (I2.4U|) for p?-i D (l - C A -i D , , x), 



< E 



a /-i a -y-a, i, a 



< 1 



Cl 



<5 A - lD (^) Q / 2 J a (X r a a 1 ,x); r° 



un\- 1 D 



un\- 1 D 



< 1 



< ci((7r 1 )- d - a + MP(7r 1 y d -P)5 x -r D (xT/ 2 F y (r^ lD < l) 

< c 1 ((7r 1 )- d ~ a + M^7r 1 y d -f 5 )5 x - lD (xr/ 2 . 



This completes the proof for (|2.41j) and hence the theorem. 



□ 



Theorem 2.8 Assume that M > and that D is an open set that satisfies the uniform ex- 
terior ball condition with radius ro > 0. For every T > 0, there exists a positive constant 
C\s = Ci$(T, ro, a, (3, M) such that for every a G [0, M\, t G (0, T] and x,y G D, 

PUt, *, y) < C 1S (i a (\ A M^! j ( t -a /o A t ja(X) y) ) (2 42) 

Proof. Fix T,M > 0. By Theorem 12.71 there exists a positive constant c\ = ci(T, ro, ct, /3, M) 
such that for every t G (0, T], 

Pt-i/ aD (l,x,y) < Cl [lAJ at (x,y))5 t -i /aD (xn - (2-43) 

Thus by ([L"5]) . (Oj) and (12^31) . for every f < T, 

< dr"/" (l a ^^(rV-x,*- 1 /^)) ^v^r 1 /^) / 2 

= ci (V d/Q AtJ a (x,y) N 



M*)° /2 
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By symmetry, the above inequality holds with the roles of x and y interchanged. Using the semi- 
group property for t < T, 

p a D (t,x,y) = [ p a D (t/2,x,z)p a D (t/2,z,y)dz 

< c 3 / p (t/2,x,z)p (t/2,z,y)dz 



t 



D 



< 9 fol ^" r % ,»). 

This proves the upper bound (|2.42p by noting that 

(1 A a)(l A b) = min{l,a, b,ab} for a,b > 0. 

□ 



3 Lower bound estimate 

Lemma 3.1 For any positive constants A, k and b, there exists C19 = Cig(A, k, b, a, /3, M) > 
such that for every z G M d , A G (0, A] and a G (0, M], 



Proof. By [HI Proposition 4.9], there exists e = e(A, k, a, (5) > such that for every A G (0, A], 



Suppose 6 > e then by the parabolic Harnack principle in [8, Proposition 4.12] 

c iP B (y,K\ 1/a )( £X,y,W ' ) < p l B{y KXl/a) {b\y,w) for w G B(y,K\ 1/a /2), 
where the constant c\ = c\ (k, b, a, (3, A) > is independent of y G M. d , A G (0, A]. Thus 

¥ v ( T B( y ,K\W) >b\)= p l B{ 1/a) (b\,y,w)dw 



> 



PB(»,«Al/«)( 6A .f.« ; ) dw 
B(y,n\ 1 / a /2) yy ' 



> ci / pL VK \v°'/2)(. EX >v> w ) dw - Cl / 2 - C 3 - 1 ) 

JB(y,K,X l / a /2) ^' ' ' 
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For the general case, by (jl.3p and (|3.1|) . 



\y-z\<K\ 1 / a 

> P (^(0,«AX/«) > &A 

^Bfo KX i/ a) (bX,0,w)dw 

B(0,KX 1 / a ) K ' ' 

p (a a -i 3 bX,0, z)dz 



= P (r 1 (B{0,K\ 1/a a ^))>a^b\\ 

> P (V(5(0,KA 1/a M~^)) > M^&a) > c 2 (A,K,6,a,/3,M) > 0. 
This proves the lemma. □ 



Recall that ip a is defined in (EHHl) . 

Proposition 3.2 Suppose that M,T > and (t,x,y) G (0,T] x D x D u»t/i S D (x) > t x l a > 
2\x — y\tp a (\x — y|) _1// ^ +a \ Then there exists a positive constant C 2 q = C 2 o(M, a, (3, T) such that 
for all a G (0, M) 

p a D (t,x,y) > C 2Q t- d ' a . (3.2) 



Proof. Let t G (0,T] and x,y £ D with > t 1 /" > 2\x - y\^ a (\x - y])- 1 /^) . By 

the parabolic Harnack principle in [8, Proposition 4.12] and the scaling property, there exists 
ci = ci(M, a, (3,T) > such that for all a G (0, M], 

p a D (t/2, x, w) < ci p a D (t, x, y) for u; G 2t 1/o /3). 

This together with Lemma 13. II yields that 

p a n (t,x,y) > — : — . , — — - / p a n (t/2,x,w)dw 
° y 1 ~ c x \B{x,tV<*/2)\ Jb^«/2) Dy ' ' 

> c 2 t- d/a [ p a B(xtl/a/2 Jt/2,x,w)dw 

JB(x,t 1 / a /2) y ' 11 

= c 2 t- d / a F x {r^ tl/a/2) >t/2) > c 3 r d /«, 
where a = a(T, a, (3, M) > for i = 2, 3. □ 



Lemma 3.3 Suppose that M, T > ; D is an open subset ofR d and (t,x,y) e(0,T]xDxD with 
min{<5i)(a;), Sd(u)} > t 1 ^" and \x — y\ a > 2~ a tifj a (\x — y|) a /( d+Q ). Then there exists a constant 
C 2 i = C 2 i(a, (3, T, M) > such that for a G (0, M] 

P x (X? D G B(y, 2- l t 1 / a )) > C 21 t d / a+l J a (x,y). 
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Proof. For t G (0,T], it follows from Lemma 13.11 that, starting at z G B(y, 4~ l t l / a ), with 
probability at least c\ = c\(a, (3, T, M) > 0, for any a G (0, M], the process X a does not move more 
than Q^t 1 ^ by time t. Thus, it suffices to show that there exists a constant c 2 = 02(0, (3, T, M) > 
such that 

'x a > D hits the ball B(y, 4T l t l/a ) by time i] > c 2 t d/a+1 J a (x,y) (3.3) 



for all a G (0,M], t G (0,T] and |s - y\ a > 2' a ti) a {\x - y \) a /( d + a ). 

Let B x := B(x, 6 _1 t 1//Q ), £ y := fl(y, G^ 1 /") and r£ := r^. It follows from Lemma ED that 
there exists c 3 = c 3 (a, /?, T, M) > such that for a G (0, M] and i G (0, T], 



E x [t A t£] > tF x (r a x >t)> c 3 i. 

By the Levy system in (|2,25p . 

P^ (x a ' D hits the ball B(y, 4~ 1 t 1 / Q ) by time i 
> P x (X? AT a G B(y, 4~V /a ) and t A r" is a jumping time 



> E 7 



J a (X a s ,u)duds 



Note that 



|x - y\ > 2-H 1 / a ^ a (\x - y\) 1/(d+a) > 2-H l / a . 
Moreover, if s < t® and u G B y , 

\X% - u\ <\x - y\ + \x - X^\ + \y-u\< 2\x - y\. 

Thus from fl33|) we get that for any a G (0, M) and t G (0, T], 

F x (x a > D hits the ball B(y, 4~ 1 i 1/a ) by time tj 

> E x [tAT%] I j a (2\x-y\)du 

J By 

> c 4 t\B y \f(2\x-y\) > c 5 t d / a+1 f(2\x-y\) > Cb2 ' d - a t d ' a+ ^ a 



J [\x-y 



for some positive constants c, = Ci(a, (3, T, M), i = 4,5. Here in the second inequality, 
□ 



(3.4) 



(3.5) 



is used. 



Proposition 3.4 Suppose that T > 0, M > 0, D is an open subset of M. d and (t, x, y) G (0, T] x 
D x D with mm{5 D (x), 5 D (y)} > {t/2f/ a and \x - y\ a > 2~ a - 1 t^ a (\x - y \) a /( d + a ) . Then there 
exists a constant C 22 = C 22 (a, P, T, M) > such that for any a G (0, M], 

p a D (t,x,y) > C 22 tJ a (x,y). (3.6) 
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Proof. By the semigroup property, Proposition 13.21 and Lemma 13.31 there exist positive constants 
ci = ci(a, (3, T, M) and C2 = c 2 {a, f3, T, M) such that for any t G (0, T] and a G (0, M] 

p a D {t,x,y) = [ p a D (t/2,x,z)p a D (t/2,z,y)dz 

> f pUt/2,x,z)p a D (t/2,z,y)dz 

JB(v,2-l(t/2)V<*) 

> Cl r d ' a T? x (x$ G B(y, 2" 1 (t/2) 1/a )) 

> c 2 tJ a {x,y). 

□ 



In the rest of this section, we assume that D is an open set in M. d satisfying the uniform interior 
ball condition with radius tq > in the following sense: For every x G D with 5o(x) < tq, there is 
z x G 3D so that \x — z x \ = 5d(x) and -B(xo> r o) C D for xq := z x + tq(x — z x )/\x — z x \. Clearly, a 
(uniform) C ' open set satisfies the uniform interior ball condition. 

The goal of this section is to prove the following lower bound for the heat kernel p%(t, x,y). 

Theorem 3.5 For any M > and T > 0, there exists positive constant C23 = (^23(0;, (3, T, M, ro) 
such that for all a G (0, M] and (t,x,y) G (0,T] x D x D, 



P%(t, x, y) > C 23 [l A j [l A j (t^ A t J«(x, y) ) . 

To prove this result, we will first prove a lower bound estimates on the Green function of X a,u 

/•oo 

G<fj{x,y) ■= / Pu(t,x,y)dt 



when U is a bounded C ' open set. The tool we use to establish the Green function lower bound is 
a subordinate killed a-stable process in U. We first introduce this subordinate killed process first. 

Assume that U is a bounded C 1 ' 1 open set in M. d and Ri the radius in the uniform interior and 
exterior ball conditions. Then it follows from [3j Theorem 1.1] that the killed a-stable process X u 
on U has a density pjj(t,x,y) satisfying the following condition: for any T > there exist positive 
constants C2 > C\ depending only on a, T, Ri and d such that for any (t, x, y) G (0, T] x U x U, 

ro(t , I , 9)aci ( 1A ^)( 1A M^)( t ^ A _i^), ( , 7) 

Let {T" : i > 0} be a subordinator, independent of X a , with Laplace exponent 

Then the process {Z^' U : t > 0} defined by Z^' U = Xj, a is called a subordinate killed stable process 
in U . Since 4> a is a complete Bernstein function, the subordinate T a has a decreasing potential 
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density u a (x). In fact u a (x) is completely monotone. (See [2TJ[25] f° r the details.) Then it follows 
from [25] that the Green function Ry(x,y) of Z a,u is given by 



R a u(x,y) 



Pu (t,x,y)u a (t)dt. 



(3.9) 



It follows from [26] that the Green function of X a,u and the Green function of Z a ' u satisfy 
the following relation: 

R a u(x,y) <G a u(x,y) (x,y)eUxU. (3.10) 

So we can get a lower bound on Gy{x, y) be establishing a lower bound on Rfj{x, y). The following 
result gives sharp two-sided estimates on R^(x,y) and the idea of the proof is similar to that of 



Theorem 3.6 Suppose that M > and U is a bounded C 1,1 open set in U. There exist positive 
constant C25 > C24 depending only on (a, /3,d,R\,M, diam(U)) such that for all a G (0, M\, 



\ \x-y\ a / 



x-y 



a—d 



log 1 + 



Su(x) a/2 Su(y) a/2 



(^(x)^(y))' r -" /2 A ^^)|^^ a/2 



and 



R a u{x,y)<C 2 A 



(5[/(z)<5[/-(y)) 



(a-l)/2 A ^(^)" /2 ^fa) a/2 



u>/ien d > a, 
when d = 1 = a, 
when d = 1 < a, 

when d > a, 
when d = 1 = a, 
when d = 1 < a. 



Proof. Since the drift coefficient of T a is 1, we know that u a (t) < 1 for all t > 0. Now the upper 
bound on Rfj follows immediately from (|3.9[) and Corollary 1.2]. Thus we only need to prove 
the lower bound. 

By using a scaling argument, one can easily check that 

u a (t) =u 1 (a^t) t>0. (3.11) 

Let T = diam(C7). Since «*(£) is a completely monotone function with it 1 (0+) = 1, by (|3,lip . 



u a (t) > u X (M^=?T) for every t G (0,T] and a G (0,M\. 
Using (I3TT2|) . (I3T91) and [3, (4.2)] we get that 



(3.12) 



1 A 



,a/2 



> 



v Vt 

u l (M~^T) 
\x — y\ d ~ a 



1 A 



$u(y) 



a/2 



-d/a 



A 



\x-y 



d+a 



u a {t)dt 



Ua 2 A U 



-3 



1A Vjggfc)^ ( 1A Vu6u(y) a/ " 



\x-y 



a/2 



\x-y 



a/2 



du. (3.13) 



Now we can follow the proof of [3] Corollary 1.2] to get the desired lower bound. In fact, when 
d > a, the desired lower bound follows from f)3. 13|) and O (4.3) and (4.7)]. Let 

duixY/HuiyTl 2 



uq :-- 



\x - y\ 
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When d = a = 1, by (13T3D and [H (4.3) and (4.9)], 



R a u{x,y) > u l (M«-PT) / Pu (t,x,y)dt 



r 



o 



> ci ) 1 A I 1 A I + C1 log(„„ V 1) + C1 « C (!/«) A 1 - J£ZlE 



\x - y\ a/2 I \ \x - y\ a l 2 



> c 2 (l Am ) + c 2 log(n V 1) +c 2 u ( (l/«o) A - 



> c 3 (l A it ) + c 3 log(u VI) > C4 log 1 + 



_ |x - 
T 

^(^) a/2 ^(y) Q/2 



Lastly, in the case d = 1 < a < 2. By (|3.13p . [31 (4.3) and (4.7)] and the first display in part (iii) 
of the proof of O Corollary 1.2], we have 

/>oo 

Ru&y) > u\M^T) I Pu {t,x,y)dt 
1 



T 

1 

> CB] 

\x-yr " 



^ (((«o V l) 1 -^) _ i) + C5U0 \^ U0 v l)-V- - (!f-Jd! 

^ 1 / i-(i/<*A 



C6 (M^OM?/)) " A : : 

V \x-y\ 



□ 



By integrating the lower bound in Theorem 13.61 with respect to y and applying (|3.10p . we obtain 
the following lower bound on E x [tj}] 

Corollary 3.7 Suppose that M > and U is a bounded C ' open set in U. Then there exists a 
constant C 2 6 = C 2 6(a, /?, d, M, R\ , diam(U)) > such that for every a S (0, M] and x £ U, 

E x [t%] > C 26 5u(xr/ 2 . 

We will first establish Theorem 13.51 for small T, that is, we will first assume that 

t<T :=(^Y. (3.14) 
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By integrating (|1.3p with respect to t and y, we have that for every open set U, A > and x E U, 

W] = / G a u (x,z)dz = \ a [ Gf:-* )/0 (\- l x,y)dy = \ a E x -i x \r aXia ~ 0)/ " (X-'U)} • (3-15) 
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Lemma 3.8 Suppose that M > 0, « G (0,1) and f/iaf G (0,T ] x D with 5 D (x) < 3t 1/Q < 

ro/4. Let z x G clD 6e smc/i i/iai \z x — x\ = Sd(x) and define n(z x ) := (x — z x )/\x — z x \. Put 
x i = z x + 3t l < a xi(z x ) and B = B(xi,3t 1 ^ a ). Suppose that xq is a point on the line segment 
connecting z x and z x + Qt l / a n(z x ) such that B(xo,2Kt 1 ^ a ) C B \ {x}. Then for any b > 0, there 
exists a constant C 2 j = C 2 -/(k, a, fl, to, b, M) > such that for all a G (0, M] 

X a h f G B{x^Kt l ' a )) > C 27 t~ l l 2 5 D {xr' 2 . (3.16) 

Proof. Let < K\ < k and assume first that 2~^K\t v l a < 5d(x) < 3t l l a . Repeating the proof 
of Lemma [3T3l we get that, in this case, there exists a constant c\ = ci(a, /3, «i, M, ro,b) > such 
that for all a G (0, M] 

F x [x^ D G B(x , Kit 1/Q )) > Cl t d / a+1 J a (x,x ) > Cl A(d,-a)t d/a+1 \x - x \- d - a 

for all t < T . Using the fact that \x - x \ G [2nt l / a , 6t 1/a ] we get that for all a G (0,M], 

P,. fe D G B(x , M 1/Q )) > c 2 > (3.17) 



for some constant c 2 = C2(a, /3, «i, M, ro, 6). By taking ki = k, this shows that ()3.16p holds for all 
b > in the case when 2" 4 Kit 1 / Q < 8 D {x) < 3t l / a . 

So it suffices to consider the case that 5 D (x) < 2" 4 Kt 1 /«. We now show that there is some 
bo > 1 so that (|3.16p holds for every b > bo and 5d(x) < 2~ i Kt 1 / a . For simplicity, we assume 
without loss of generality that xq = and let B := B(0, nt 1 ^). Let x 2 = z x + A~ l K\i{z x )t l / a and 
B 2 ■= B(x 2 ,4T 1 Kt 1 / a ). Observe that since B(0, 2«i 1 /°) c S\ {x}, 

K /2t 1/a < \y - z\ < Qt 1/a for y £ B 2 and z G 5(0, «;t 1/Q ). (3.18) 

By the strong Markov property of X a at the first exit time r^ 2 from B 2 and Lemma |3.1| 

F x (x^eB^Kt 1 /*)) 

> F x (t% 2 < bt, X«„ G B(0, 2~ 1 Kt 1 l a ) and \X a s - X^ \ < k/2 for s G [r| 2 , r£ 2 + hi 1 / ]) 

> c 3 (r£ 2 < bt and G B(0, 2" 1 ^ 1/a )) . (3.19) 

It follows from f)3. 15j) and Corollary 13.71 that 

F X (X% eBiO^- 1 ^)) = [ [ G a B2 (x,y)J a (y,z)dydz 

'-IE, [rl 2 



> 



> 



c 4 E x/tl/ . O ( r;/,,, V 1-', 

M 



> (3-20) 
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for some positive constants 04,05 depending only on a,f3,ro,K and M. Note that, by (jl.3p 

/ p a B{ lKtl/a) (bt,x,z)dz= pf^jill lK) {b,r 1 l a x,w)dw. 

Since at^ a ~^^ a ^ < MT^ a ^/ a/3 ; by applying Theorem 12,81 to the right hand side of the above 
display, we get 

> bt) < b~ d ' a / B[t X2A p dw 



t X/a 

for some positive constant cq depending only on a, (3, ro, k and M. Define 



a/2 



(3.21) 



bo 

We have by ([3TT9V([3T2Tj ) that for b > b , 



2c 6 "\ 2 rf+« 

C5 



¥ x (xSt eB) > c 3 (w x ( x h 2 G B ^ 2 " 1 ^ 1/a )) - p - ( T k > bt 

* ^(c 5 /2)(^) a/2 . (3.22) 

()3.17p and (I3.22p show that (|3.16p holds for every b > bo and for every x £ D with 5d{x) < St 1 ' . 

Now we deal with the case < b < b and 5 D (x) < 2- i nt 1 / a . If S D (x) < 3(bt/b ) 1/a , we have 
from (|3.16p for the case of b = bo that 

F x (*£ G £(x Q , Kt 1 '*)) > P x . (x° (6t/6o) G B(x , K(bt/b 

( Spjx) V /2 _ c (Sn(x)^ ' 



If 3(fa/6 ) 1/a < o D (x) < 2~ 4 Ki 1 / a (in this case k > 3 • 2 4 (6/6 ) 1/a ), we get (^T6|) from (13371) by 
taking rci = (6/6q) 1//q . The proof of the lemma is now complete. □ 



Proposition 3.9 Suppose that M > and G (0,T ] x D x D with \x - y\ < t^t^Qx - 

y|) 1 /( a! + a ) ; 5 D (x) < 2t l / a and 5o{y) < ro/5. T/ien i/iere exists a constant C28 = C , 2s(a, /3, M, ro) > 
swc/i t/iat for all a G (0, M], 

p a D (t,x,y) > C 28 ^^ 1 o D (x) Q / 2 o D (y) Q / 2 . (3.23) 



Proof. Under the assumptions of the proposition, there are points z x , z y G dD and xo, yo G D such 
thatfe(x) = \x-z x \, 8 D {y) = \y - z y \, 8B(x , 4t l/a ) n dD = {z x } and dB(yo,At l l a ) f] dD = {%}. 
Observe that 

S D (x ) = S D {y ) = 4iV* an d |x - x |, |y - y | G [* 1/a , ^ a ). 
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By the semigroup property, with B := B(xq,A 1 t 1 ^ a ) and B := B(yo,4 l t l l a ) 
p a D {t,x,y) = f p a D {t/3,x,z) f p a D (t/3,z,w)p a D (t/3,w,y)dwdz 

JD J D 

> [ p a D {t/3,x,z) fp a D {t/3,z,w)p a D (t/3,w,y)dwdz 

J B J B 

> inf ^p a D {t/3,z,w) / p%(t/3,x,z)dz (p a D (t/3,w,y)dw. 

(z,w)eBxB Jb JB 

Since for z £ B and w G B, 

5d(z) > 5 D (x ) - \x -z\> t 1 ^, 8 D (w) > 5 D (y ) - \y -w\> t l ' a , 

\z-w\<\z- x \ + \x Q -x\ + \x-y\ + \y- y \ + \yo - w\ < lOt 1 ^^ - y\) 1/{d+a \ 

by combining Proposition l3.2l and Proposition l3.4l we have that there exists c\ = c\(a, f3, ro,M) > 
such that for all a G (0,M], 

inf _p%(t/3,z,w) > Cl r d / a . 

(z,w)£BxB 

Since Sd(x) < 2t x l a < tq/8 and 5d{u) < 3t 1 l a , we have by Lemma I3T81 

p a D (t,x,y) > c 2 r d l a - 1 5 D (x) a ' 2 5 D (yr/ 2 
for some positive constant C2 = 02(0, (3, M, ro). □ 

Proposition 3.10 Suppose that M > and (t,x,y) G (0,T ] x D x D with 5 D (x) < (t/2) 1 / < 
o~d{v) and \x — y\ a > t^ a (\x — y\) a /( d+a ) . Then there exists a constant C29 = ^29(0, (3, M, ro) > 
such that for all a G (0, M], 

p a D (t,x,y) > ^t^Soixr^JO&y). (3.24) 

Proof. Since 5 D (x) < (t/2) 1 / < r /16, there are z x G dD and zq G D such that 5d(x) = \x — z x \ 
and dB(z 0l 2t 1 ' a ) n dD = {z x }. Choose x in Bfo, 2i 1 / a ) an d K = K ( a ) e (0, 1) such that 

B (x , 2«tf 1 / a ) C 5 (* , (2 - 2" 2 / Q )t 1 / Q ) n 5 (z, (1 - 2- 1 - 2 / a )t 1 ' a ^ . 

Such a ball B(xq, 2Kt 1 > a ) always exists because 

2 < (2 - 2~ l ) + (1 - 2- 2 ) < (2 - 2~ 2 / a ) + (1 - 2" 1 - 2 /"). 
Since \x - y\ > t l / a tp a {\x - y\) l/(d+a \ we get that for every z G B(x , nt l / a ), 5 D (z) > (t/4) 1/o and 

\y-z\>\y-x\-\z-x\> 2~ 1 (i/4) 1 / c > a (|x - y\) 1/{d+a) . 
On the other hand, for every z G B(xo, Kt 1 /"), 

\z-y\ < \z-x\ + \x-y\ < (l-2~ 1 ' 2/a )t 1/a + \x-y\ < 2\x - y\. 
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Thus by the semigroup property and Propositions 13.21 and 13.41 there exist positive constants c±,C2 
and C3 depending only on (a, (3, tq, M) such that for all a G (0, M], 



p a D {t,x,y) = f p a D (t/2,x,z)p a D (t/2,z,y)dz 
Jd 



> / p a D (t/2,x,z)p a D (t/2,z,y)dz 

JB{x Q ,Kt 1/a ) 



> at / p a D (t/2,x,z)J a (z,y)dz 

JB(x ,Kt 1/a ) 

> c 2 tf(2\x-y\) [ p a D (t/2,x,z)dz 

JB(x ,Kt 1/a ) 

> c 3 tf(\x - y\) F x (X a /2 D G B(x , nt l ' c 



t/2 

Applying Lemma 13.81 we arrive at the conclusion of the proposition. □ 

Proposition 3.11 Suppose that M > and (t,x,y) G (0,To] x D x D with 

m<xx{5 D (x), 8 D {y)} < (t/2) 1 /" < \x - y\r{\x - y\)~ 1 /^ . 
Then there exists a constant C30 = Cso(a, (3, M, ro) > such that for all a G (0, M], 

p a D (t,x,y) > C- i0 5 D (x) a ' 2 5 D (y) a ' 2 J a {x,y). (3.25) 

Proof. As in the first paragraph of the proof of Proposition 13.91 set z x G dD and xq G D so 
that \x - z x \ = 5 D (x) and dB(xo,?>t l l a ) n dD = {z x }. Let k := 1 - 2-~ i / a . Note that for every 
z G B(xq, Kt l / a ), we have 

4t 1/a > 5 D (z) > 2(t/2) 1 / Q . 
If \y-z\< t l / a ijj a {\y - zl) 1 /^^, we can apply Proposition 13.91 and the assumption 



(t/2) 1 ^ <\X- y\lp a {\x - y\)-V(d+a) 



to get that 



p a (t/2,z,y)>c 1 t 1 / 2 5 D (y)^ 2 J a (x,y). 
If \y-z\> t l l a ^ a {\y - z]) 1 /^^, we can apply Proposition 13.101 to get that 

p a (t/2,z,y)> c 2 t 1 / 2 5 D (yT/ 2 J a (y,z). 

For z G Bixo,^ 1 ^), we have 

\z — y\ < \x — y\ + |xo — x\ + |xo — z\ < \x — y\ + 4t 1//c 

< |z - y| + 2 2+1 / Q (t/2) 1 / a V a (k - yl) 1 /^) 

< (l + 2 2+1 / a )|x-y|. 

Thus if |y - z| > ^/"^(ly _ z |)V(d+a) ) we have 

p a (i/2 )2) y)>c 3 t 1 / 2 <5 i3 (2/) a / 2 J a (x ! y). 
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Consequently we have for all z S B(xq, nt 1 ^") 

p a (t/2,z,y) > c 4 t 1 ' 2 5 D (y) a / 2 J a (x,y). 
Hence by the semigroup property we get 

p a D (t,x,y) = [ p a D (t/2,x,z)p a D (t/2,z,y)dz 

JD 



> / p a D (t/2,x,z)p a D (t/2,z,y)dz 
JB(x Q ,Kt 1/a ) 

> c 4 f p a D (t/2,x,z)t^ 2 S D (y) a / 2 J a (x,y)dz 

J B{x ,Kt 1/a ) 

= c A t 1 / 2 6 D (y) a / 2 J a (x,y) [ p a D (t/2, x, z)dz 

= c 4 t 1 / 2 5 D (y) a / 2 J a (x,y)P x (xjf G B(x Q , nt 1 / 01 )) 
= c 5 5 D (xr/ 2 5 D (yr/ 2 J a (x,y). 
We arrive at the conclusion of the proposition. □ 

Proof of Theorem 13.51 In this proof, for two non-negative functions / and g, the notation 
/ x g means that there are positive constants c% and C2 depending only on M, d, a and (5 so that 
c\g{x) < f(x) < C2g(x) in the common domain of definition for / and g. 
We first assume that t < Tq. 

1. We first consider the case \x — y\i[) a (\x — y|) _1 /( d + Q ) < t 1 /". We claim that in this case 

fB(ll , 1 , ) >^( 1A ^!)j 1 ^j. (3.20) 

This will be proved by considering the following two possibilities. 

(a) max{5£)(x), Sd(u), \% — y\4> a (\x — y\)^ l ^ d+a ^} < t 1 /": Proposition 13.91 and symmetry 
yield (1336!) 

(b) m&x{5 D (x),5 D {y)} > t x l a > \x - y\^ a (\x - : 

lfm&x{5 D {x),8 D (y)} > t x l a > 2\x-y\if; a (\x-y\)- 1 / ( - d+a \ (l3T25j) follows from Proposition 



If mm{5 D (x),5 D {y)} >t l / a and \x - y\ < t x l a < 2\x - y\^ a (\x - y|)~V(<*+a), 

m\x- y \) ^ d/a ( Sp(xr/ 2 \ ( 8 D (yr/ 2 \ 
\x-y\ d+a \ V~t J \ VI J' 

If max{(5£> (x),S D (y)} > mm{5 D {x), 5 D (y)} < t l / a and \x - y\ip a (\x - y^ 1 /^) < 
t 1 ^ < 2\x - y\ip a {\x - y\y l/{d+a) , 

Thus by combining Proposition 13.41 and Proposition 13.91 we get (|3.26p for the case of 
max{5 D (x),6 D (y)} > t 1 ^ and \x - y\ < t l / a < 2\x - y\tp a (\x - y\)- 1 K d+a h 
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2. Now we consider the case \x — y\ip a (\x — y\) 1 /i d+a ) > t 1 /" and claim that 

PD(t , x , y) >J 1A Sj^\ ( 1A ^l\ tf(}x . yn . (3.27) 

(a) mm{5 D (x),5 D (y)} < {t/2) l / a and \x - y\ip a (\x - > t x l a : By symmetry we 
can assume 8r>{x) < (t/2) 1 /". Thus combining Propositions 13.101 and 13.111 we have 
(|3,27[) for this case. 

(b) mm{5 D (x),8 D {y)} > (t/2) 1 / a and |x - y|V> a (|x - y|)~ 1/(a!+a) > t l / a . In this case, clearly 

Thus Proposition 13.41 yields (|3.27p . 
We have arrived at the conclusion of Theorem 13.51 for t < Tq. 



Assume T = 2T . Recall that T = (r /16) a . For (t,x,y) e (T ,2T ] x D x D, let x ,y € 
be such that max{|x — xo\, \y — yo\} < ro and min{<5£)(xo), Sd(uo)} > ro/2. Note that, since 

< cii a (2r), for all r > 0, (3.28) 

if | a; — y| > 4ro, then — y\ < \x — y\ — 2r$ < \x$ — 2/0 1 < \ x — y\ + ^ r o < fl^ — 2/|> and so 
c^ 1 J a (xo, yo) < J a (x,y) < C2J a {xo,yo) for some constant C2 = C2(M) > 1. Thus by considering 
the cases \x — y\ > 4ro and |x — y| < 4ro, we have 

{t/2yd/a A tja (^yo) > C3 / r<i /a A (t ja(x5 (3 29) 

Similarly, there is a positive constant C2 such that 

(V3)-^A^) > c 4 ((t/(12))-V« A ^^1) , z e D, 

W)- d/a A Ua{ ^ y) > c 4 ((V(12))-^ A V w £ D. (3.30) 

By (|3.30p and the lower bound estimate in Theorem 13.51 for p^ on (0, To] x D x D, we have 



p a D {t,x,y) = I p a D (t/3,x,z)p a D (t/3,z,w)p a D (t/3,w,y)dzdw 



V~t J V JJdxd\\i2J 12 y ^ y^/3 

p&(tA*,«>)( (jgj A — 12 ) I vP~J 
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for some positive constants Q,i = 3,4. Let D\ := {z G D : 5d(z) > ro/4}. Clearly, xo,yo £ D\ and 



min{5 Dl (xa),5 Dl (yo)} > r /4 = 4(T ) 1 / a > 4(i/2) 1 A*. (3.31) 
By (fl~4"|) and (f3T29j) . we have 

Dxfl 



12; 12 j y v /t/3 y 



>c 8 / p a (i/(12),x ,z)pf )l (t/3,z, 1( ;> a (t/(12), U ;,yo)d^ 



>c 8 / ^ 1 (t/(12),x ,z)p^ 1 (i/3,z,u;)p^ i (t/(12), U ;,y )^ ( i U ; 

= c 6 p a Dl (t/2,x Q ,y ) > c 9 r(V2)- rf / Q A ^"'^ "j > Cl0 A (tJ a (x,y))) 



for some positive constants Ci,i = 7, • • • ,10. Here both Propositions 13.21 and 13.41 are used in the 
third inequality in view of (|3.3ip . By repeating the argument above, we have proved Theorem 13.51 
□ 

Proof of Theorem 11.11 Theorems 12.81 and 13.51 give Theorem I l.lf i). By [15} 118]. for any bounded 
open set D in R rf , X a ' D is intrinsically ultracontractive. Since the function tp a (\x — y\) is bounded 
above and below by a positive constant if D is bounded, using the intrinsic ultracontractivity of 
X a on bounded open set and the continuity of eigenvalues proved in |10j . the proof of Theorem ll.il 
(ii) is almost identical to the one of [5j Theorem 1.1(h)] . We omit the details. □ 



Proof of Corollary 11.21 The lower bound estimate in (jl.5p follows from (|3,10p and Theorem! 

Since the function ip a {\x—y\) is bounded above and below by a positive constant if D is bounded, 
by integrating the two-sided heat kernel estimates in Theorem 11.11 with respect to t, the proof of 
the upper bound estimate in (jl.5p is identical to the one of Corollary 1.2] so we omit its details 
here. □ 



Theorem 3.12 (Uniform boundary Harnack principle) Suppose M,R G (0, oo) and that D 
is an open set in W 1 , z G dD, r G (0, R) and that B(A,nr) C D n B(z,r). There exists C31 = 
Czi(d, a, P, k, M, R) > 1 such that for every a G (0, M], and any functions u, v > on M. d , positive 
regular harmonic for X a in D n B(z, 2r) and vanishing on D c n B(z, 2r), we have 



„ 1 u(A) u(x) „ u(A) 
-y(j4j v[x) u(A) 
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Proof. Note that by the approximate scaling property in (jl.3p . we have for every r > 0. 

G a Bi0 , r) (x,y) = r a - d G a ^ )/9 (x/r,y/r). (3.32) 

Thus applying [3, Corollary 1.2] and our Corollary 11.21 to (|3.32p . we have that for every R, M > 0, 
there exists c = c(a, (3, R, M) > such that, for every a £ (0, M] and < r < R 

c~ 1 G ! B(xo,r)(a; ) 2/) < G B(x ,r)( x >y) - cG B(x ,r)(x,y), Vx, y G B(xo,r). (3.33) 
Using (I3.33h . we can get uniform estimates on the Poisson kernel 

K B(x ,r)( x i z ) ■= / G a BM (x,y)J a (y,z)dy 

J B(xo,r) 

of B(xo,r) with respect to X a for r G (0, i?]. In particular, for r < |z — xq\ < 2R, Kg, r Jx,z) is 
comparable to K B r X0>r \(x, z), the Poisson kernel of B(xq, r) with respect to X for r £ (0, i?]. Then 
using the uniform estimates on K b , Xq r \(x, z) and (|3.33[) we can easily see that [271 Lemma 3.3] can 
be proved in the same way. Using the uniform estimates on the Poisson kernel of B(xq,v), f|3.28[) 
and (|3.33p we can adapt the argument in [T|ll91[27j to get our uniform boundary Harnack principle. 
We omit the details. □ 

Proof of Theorem 11.31 First we observe that Harnack inequality holds for the process X := X 1 
by [21]. That is, there exists a constant c\ = c\(a, /3,M) > such that for any r £ (0, M^^ a ~^], 
Xo £ M. d and any function v > harmonic in B(xo,r) with respect to X, we have 

v (x) < c\v{y) for all x, y £ B(xq, r/2). (3.34) 

Note that for any a £ (0,M], X a has the same distribution as {XXx~a t ,t > 0}, where A = 
a P/(P- a ) > j\^f/3/(/3~ Q ). Consequently, if u is harmonic in B(xo,r) with respect to X a , where r £ 
(0, 1], then v(x) := u(Xx) is harmonic in ^(A^xo, A _1 r) with respect to X and A~ x r < M^/^~ a \ 
So by (13341 

u(Xx) = v(x) < c\v{y) = c\u(\y) for all x,y £ B(X^ 1 xq, A _1 r/2). 

That is, 

u(x) < c\u(y) for all x, y £ B(xo, r/2). (3.35) 

In other words, uniform Harnack inequality holds (for every r < 1) for the family of processes 
{X a ,a£(0,M}}. 

Since D is C 1,1 open set, there exists tq < Rq such that the following holds: for every Q £ dD 
and r < ro there is a ball B = B(zq,r) of radius r such that B C M. d \ D and dB n dD = {Q}. 
In addition, it follows [231 Lemma 2.2] that, for each Q £ dD, we can choose a constant c% = 
C2(d, A) £ (0, 1/8] and a bounded C ' open set Uq with uniform characteristics (J?*, A*) depending 
on (.Ro, A) such that B(Q, c 2 r ) n D C Uq C B(Q, r ) n D and 

<^£»(y) = ^C/ Q (y) for every y G B(Q, c 2 r ) n D. (3.36) 

Assume a £ [0,M], r £ (0, C2?"o], Q G dD and u is nonnegative function in W 1 harmonic in 
D fl B(Q,r) with respect to X a and vanishes continuously on -D c n B(Q,r). Let := ZQ r °. By 
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the boundary Harnack principle (Theorem 13. 12j) . there exists a constant C3 = 03(0;, /3, a, Rq, A, M) 
such that 

u(x) 

-j4 < cst^t r for every x, y € B(Q, r/8) D £>. 

Now applying Corollary 11,21 to Gfj (x, zq) and Gfj (x, zq), then using (|3.36p . we conclude that 



S < «^ = CiC^ f - -ery x, y G c 2 r) n I? (3.37) 

lyj 5 u Q \y) 5 d \y) 

for some C4 = 04(0, /3, a, i?c>; A, M) > 0. 

Now Theorem 11.31 follows from the uniform Harnack principle in (I3.35P , (|3.37p and a standard 
chain argument. □ 
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